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We study the diffusion on an annealed disordered lattice with a local dynamical reorganization of 
bonds. We show that the typical rearrangement time depends on the renewal rate like f,- ~ r" with 
a 7^ 1. This implies that the crossover time to normal diffusion in a slow rearrangement regime 
shows a critical behavior at the percolation threshold. New scaling relations for the dependence 
of the diffusion coefficient on the renewal rate are obtained. The derivation of scaling exponents 
confirms the crucial role of singly connected bonds in transport properties. These results are checked 
by numerical simulations in two and three dimensions. 

PACS numbers; 5.70.Jk, 72.60.+g, 82.70.Kj 
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"Dynamic percolation" ("stirred percolation") 
was introduced as a model of transport in environments 
that evolve in time, e.g. microemulsions or polymers 
(for further applications see ||3|]). The simplest version 
of the model is defined on a d-dimensional regular lat- 
tice. Each pair of nearest neighbor sites is connected by 
a bond, which can be either conducting or insulating. We 
note p the proportion of conducting bonds. Time evolu- 
tion of the environment is achieved by a reorganization 
of bonds, defined below. Diffusion of a tracer particle 
in such a network is conveniently described by the ant- 
in-the-labyrinth paradigm B. Two basic algorithms are 
available. The "blind" ant chooses its direction randomly 
at each time step and moves only if the corresponding 
bond is conducting. The "Myopic" ant chooses among 
the conducting bonds. Both algorithms lead to the same 
scaling behavior of the diffusion coefficient. Two qual- 
itatively different dynamic percolation models appeared 
in the literature. The global reorganization model is the 
simplest. After some renewal time T^, the assignment 
of conducting bonds is updated throughout the lattice. 
The behavior of this model is well understood H , as it 
is closely related to the ordinary percolation. If \r^)rp is 
the mean square distance traveled on the quenched lat- 
tice during the time Tr, the diffusion coefficient on the 
stirred lattice will he D = {r'^)rp /2dTr. The case of lo- 
cal reorganization, which is studied in this article, is more 
realistic, because the evolution of the network is contin- 
uous. The state of a bond evolves through a Poissonian 
process with a characteristic time r. At each iteration a 
conducting bond is cut with a probability l/(pr), and a 
randomly chosen non conducting bond becomes conduct- 
ing, to insure that the proportion p of conducting bonds 
is conserved. No exact result is available for the depen- 
dence of the diffusion coefficient D on p and r, except in 
some particular one dimensional situations [p|. Approx- 
imative solutions of the problem in any dimension can 
be obtained by means of a time-dependent version of the 
effective-medium approximation developed in M. 

Here, we study the scaling of the diffusion coefficient 
D in the vicinity of the percolation threshold pc of the 
quenched network. Several different scaling formulas for 



D{p — Pc, t) were proposed in the literature. They were 
derived for models with slightly different local evolution 
rules, but the details in the local rules are not relevant 
for the critical behavior around pc Q . As discussed be- 
low, our simulation results do not support current pre- 
dictions. We derive a new scaling formula for the diffu- 
sion coefficient, which we confirm by extensive numerical 
simulations. The behavior at the percolation threshold 
is studied first, before we treat the general case of the 
behavior around Pc- 

The mean square displacement in the vicinity of pc on 
a quenched percolation network is given by M] 
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where d'^ is the anomalous-diffusion exponent, d'^ 
{2iy + fi-P)/{iy- (3/2) and 



/(^) 



const. 
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At early times, anomalous diffusion is observed. The 
crossover to a normal diffusion (if p > pc) or to a lo- 
calization regime (if p < pc) appears at a time of the 



order of tr 



b-p< 



|/3-2iy-/i 



which is the only relevant 



timescale of the problem. In the case of dynamically dis- 
ordered lattices, another timescale, related to the cluster 
rearrangement process, has to be taken into account. We 
define this typical "rearrangement" time tr as crossover 
time from anomalous to normal diffusion at the percola- 
tion threshold. It is only a function of the evolution rate 
T, and we assume a dependence in the form tr ~ r". The 
mean square displacement in the presence of dynamical 
disorder is thus described by a scaling formula depending 
on two parameters, t/tc and t/tr [l|: 



(i?2) ^ t^/<,g [(p-p,)tl/(2-+A'-/5);i/r" 



(2) 



At the percolation threshold, tc diverges, and the preced- 
ing expression reads 



{R')=t'^^'^xit/r'^): 



(3) 
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FIG. 1. Scaling function x (Eq- (tBl)) in two dimensions. 
Results obtained for r = 1.62 x 10*^ (solid line), r = 4.05 x 10^ 
(dotted line), r = 1.8 x 10^ (dashed line), r = 4.5 x 10'* (long 
dashed line), r = 5.0 x 10"^ (dot dashed line) with the blind 
ant algorithm. Asymptotic behavior y ~ a;^"^'''™ (bold long 
dashed line). 



where xiv) ^ const, as y — + and xiv) = Dy^^'^''^^" as 
2/ — > oo. The diffusion coefficient D is obtained in the 
hmit t — ^ cx) by 



D 



-afi/(2v+ii-l3) 



(4) 



Eq. (0) contains an unknown parameter a. Several val- 
ues of a were proposed in the literature. In 0, the 
problem was mapped on the continuous random walk, 
and the lower and upper bounds for a were predicted. 
In [y a = 1 is considered. The only justification for 
this value is the assumption that the global and local re- 
arrangement models have the same behavior. We have 
performed Monte Carlo simulations to evaluate a numer- 
ically. The diffusion coefficient can only be measured for 
small values of t, where the crossover time tr is small. In 
order to explore a broader range of values, we have de- 
termined a from the finite size scaling relation (pi). We 
measured (i?^) for r between 5 x 10^ and 1.62 x 10® in 
two dimensions, and between 7.8 x lO'^ and 5.12 x 10^ in 
three dimensions. In two dimensions the best data col- 
lapse with parameter a is obtained for a = 0.80 ± 0.02 
(Fig. |l]), in three dimensions for a = 0.79 ± 0.03 (Fig. 
y) . Identical results were obtained with both the myopic 
and the blind ant algorithms. 

As a matter of fact, the value of a can be evaluated 
as a function of known critical exponents using simple 
assumptions about the geometry of clusters. Clusters 
are composed of well connected blobs, interconnected by 
singly connected bonds ( "red bonds" ) [|| . If a red bond is 
cut, the cluster breaks into two parts. We argue that the 
crossover time corresponds to a removal (or addition) of a 
red bond in the region visited by the tracer particle. The 
red bonds are the only possible paths where a particle 
can escape from a blob, hence they control the diffusion. 
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FIG. 2. Scaling function x (El- (h)) in three dimensions. 
Results obtained for r = 5.12 x 10^ (solid line), r = 1.28 x 10^ 
(dotted line), r = 3.12 x 10* (dashed line), r = 7.8 x 10^ 
(long dashed line) with the myopic ant algorithm. Asymptotic 



behavior y ^ x 



l-2/d'^ 



(bold long dashed line) 



For t < tr a. particle visits on average a hypersphere of 
a diameter R ~ i^/'^™. The "network" of red bonds is 
fractal, and their number inside the hypersphere grows 
as A^rb ^ i?^/" 1^. The crossover corresponds to the 
average time for the first of A^rb red bonds to be cut. 
Hence tr - r/A^^b -^ r''™/('i™+i/i')^ 



givmg 



1/^ 



(5) 



In two dimensions, where v 
/i = 1.303 lO] we obtain a = 



= 4/3, p = 5/36 Ig and 
0.802. In three dimensions 
a = 0.81 ± 0.06 is obtained, using v = 0.88 ± 0.02, fi = 
2.003 ± 0.047 ^ and (3 = 0.405 ± 0.025 |3|. These 
values of a are in excellent agreement with numerical 
results. Relation M) predicts that a = 1 for d > 6, so 
in this limit the local and the global reorganization rules 
lead to the same scaling. 

Knowing the value of a, the complete scaling law for 
D in the vicinity of the percolation threshold can be de- 
duced from (ph. The ratio tc/tr separates two different 
regimes. In the fast rearrangement regime {tc/tr ^ 1) a 
tracer particle does not see the finitcness of cluster sizes, 
hence the scaling of D is given by (^. In the slow re- 
arrangement regime {tc/tr ^ 1) two cases have to be 
considered. For p > pc and r — > oo, known results for 
the diffusion on the quenched network [Q should be re- 
covered, hence 



D r^\p~Pc 



(6) 



For p < Pc the situation is more complicated. At 
t ^ tc 'ti tr, the network is not yet reorganized and 
anomalous diffusion crossovers to a localization regime 
on a finite cluster exactly in the same way as for the 
quenched network. The mean square displacement is 
thus (i?^) ^ \p — Pc]'^^^"- For t > tc it grows as 
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FIG. 3. Calculation of D for different values of \p — pc\ and 
T. Results for p = 0.47 (o), p = 0.43 (□), p = 0.42 (o), 
p = 0.41 (A), p = 0.40 (v), P = 0.53 (+), p = 0.56 (x), 
p = 0.60 (*), p = 0.70 (>). Function (|) (solid line) 
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For i ^ cxD a diffusive regime is reached, and it is evident 
that D ^ 1/r in this case. Thus the scaling function g' 
behaves as g'[x,y\ ~ x~°"y^i°' for a;,j/ — > oo, where the 
coefficient a reads 

a^{l/a-l){2v + ^l-P)^l~■^. (7) 

2v 

The final expression was obtained replacing a by (|5|). 
Then a = 0.948 in two dimensions and a = 0. 77 ±0.017 
in three dimensions. The scaling relation for t > tc can 
thus be written as a function of a unique parameter 



{R')^\p-p,f~''^g" 



\p~ Pc\"-T 



(8) 



with /(y) ~ const for y —>■ 0^ f{y) ~ y for y — » oo. It 
is readily seen that the crossover time t'^ = \p — Pc\'^t 
has itself a critical behavior near pc with an exponent 
a. This fact has been already predicted in [|lj but a 
different exponent a = 1 was proposed. The scaling of D 
in the slow rearrangement regime is simply deduced from 
(i, 
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The complete scaling law for D consistent with 
and (O) reads 



(9) 
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FIG. 4. Results of calculation of D for different values of 
|p — Pel and r rescaled by (|lO[). Results for p = 0.47 (o), 
p = 0.43 (D), p = 0.42 (o), V = 0.41 (A), p = 0.40 (v), 
p = 0.53 (+), p = 0.56 (x), p = 0.60 (*), p = 0.70 (o). 
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FIG. 5. Scaling function g" (ra) in two dimensions. The 
data for p = 0.4 (solid line), p = 0.41 (dotted line), p = 0.42 
(dashed line), p = 0.43 (long dashed line), p = 0.44 (dot 
dashed line) and p = 0.46 (solid line) Asymptotic behavior 
y ^ X (bold long dashed line) 



To verify this relation, we have calculated the diffu- 
sion coefficient in two dimensions for different values of 
T and for p in the range [0.4; 0.47] and [0.53; 0.7], using 
the algorithm of the myopic ant. Results are presented 
in Figure 0. They are well rescaled by the relation ( p^ ) 
(Fig. ^). The best collapse seems to be reached for a 
slightly smaller value of a (a = 0.9) than predicted by (0) 
(a = 0.948). However, the collapse is not very sensitive 
on the precise value of a, because the slow rearrangement 
regime is not explored in our range of (r, p) . It is difficult 
to attain this regime using a simple random walk, since 
the crossover time to the diffusive behavior becomes too 
important for large values of r. This is the reason why 
we used the following algorithm to verify Eq. (0). 

We start from a given site belonging to a cluster of s 
sites. We suppose that the evolution of the network is 



quasistatic : before the network is rearranged, the par- 
ticle is thermahzed, so that the probabihty to find it on 
a given cluster site equals l/s. Thus we assign at first 
the probability l/s to each cluster site. We then ex- 
change one conducting bond with an insulating bond, 
find a new cluster distribution and thermalize the prob- 
ability distribution on each cluster. We iterate this pro- 
cedure and measure the mean square displacement. The 
Hoshen-Kopelman algorithm |l3| was used to obtain the 
distribution of clusters. To get good statistics, an aver- 
age over more than 2000 realizations was performed, so 
we were limited to networks of relatively small size (up to 
400 X 400 sites) . Since the diffusive regime is not attained 
on such a small network, we used the finite size scaling 
formula (||). We measured (R"^) for p ranging from 0.4 
to 0.46. For higher values of p, clusters are too large, and 
much bigger networks have to be used. The data collapse 
is obtained for a = 0.87 ± 0.05 (Figure |^), that is for a 
value slightly smaller than predicted by (0). The same 
effect as in the case of the data collapse oi D{p~pc, t) is 
thus encountered. The discrepancy is due to fact that we 
are already out of the critical region, so corrections to the 
exponents a and d^ should be taken into account. For 
values of p far from pc, the probability of having a large 
cluster, corresponding to a long jump, grows more slowly 
than near pc, and the growth of the diffusion coefficient 
with p is thus also slower. 

In conclusion, we have derived a new scaling law for 
the diffusion coefficient in the case of a simple model of 
stirred percolation. The dependence of the scaling ex- 
ponents on the basic exponents of the percolation the- 
ory was found. We showed that the distribution of red 
bonds controls the transport in the network. Results are 
supported by extensive numerical simulations. In the 
slow rearrangement regime for p < pc the diffusion coef- 
ficient scales a.s D ^ \p — Pel" , where s' = —2.1 in three 
dimensions. The value of the scaling exponent in mi- 
croemulsions (s' = 1.2) [[l7|-[19[ thus cannot be explained 
by this simple model, as suggested earlier 16 1. It is plau- 
sible that the difference is due to interparticle interac- 
tions present in microemulsions. They play an impor- 
tant role in the formation of clusters pG| and they might 
also influence the dynamics of the reorganization of the 
environment. 
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